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We demonstrate the strong coupling between an electron spin ensemble and a three-dimensional cavity in a 
reflection geometry. We also find that an anticrossing in the cavity /spin spectrum can be observed under 
conditions that the collective coupling strength Qc is smaller than the spin linewidth 7^ or the cavity linewidth. 
We identify a ratio of Qc to 7^ {gdls > 0.64) as a condition to observe a splitting in the cavity frequency. 
Finally, we confirm that gc scales with a/TV, where N is the number of polarized spins. 



When two oscillators couple strongly, their otherwise 
degenerate modes repel each other to develop an anti- 
crossing in the spectrum. A canonical example is the 
interaction between a single photon and a single atom, 
and this phenomenon is called the vacuum Rabi split- 
ting.^ The strong light-matter coupling has been a cen- 
tral subject in cavity quantum electrodynamics (QED), 
and has attracted increasing attention as a means to co- 
herently transfer information between a flying (photonic) 
and stationary (atomic) qubits in the field of quantum in- 
formation processing. So far, the strong coupling regime 
has been realized in various systems j^ — Even when the 
cavity-atom coupling strength gg is small, as in the case 
of natural spins, the interaction can be collectively en- 
hanced, by utilizing an ensemble of N identical two-level 
systems, to replace gg with a new parameter gc = gs^/N - 
Such a collective enhancement has been observed in a 
wide variety of systems, such as Bose-Einstein conden- 
sates^il trapped ions,^ and circuit QED with supercon- 
ductors.- 

More recently, the strong coupling between an electron 
spin ensemble and a superconducting coplanar waveguide 
resonator has been demonstrated for impurity spins in 
ruby and diamond^i^iii A rough estimation of gs in the 
spin system is given as gs ~ tyiq^ jiQUJc/'^hVc^ where mo 
is the magnetic dipole moment of the spin, ujc is the cav- 
ity frequency, and Vc is the cavity mode volume. There, 
the small Vc of the order of 10~^^ m^ helped to enhance 
the cavity-spin coupling. These experiments are moti- 
vated by the goal of using an electron spin ensemble as 
a medium for quantum memory, which will potentially 
benefit from the capability of storing multiple qubits in 
different modes of the ensemble. 

In this Letter, we study the coupling between an elec- 
tron spin ensemble and a microwave resonator with much 
larger Vc^ a three-dimensional cavity that is commonly 
used for bulk electron paramagnetic resonance (EPR), 
and demonstrate that by increasing the number of spins 



Electronic mail: eisuke. abe@materials. ox. ac.uki 



the spectrum experiences the transition from the weak to 
strong coupling regimes. Our result provides an opportu- 
nity for studying the strong coupling effect with readily 
accessible setups and wider tunabilities in the experimen- 
tal parameters. In addition, it has been recently observed 
that such a three-dimensional cavity allows a supercon- 
ducting qubit to preserve its coherence longer than the 
case of a two-dimensional counterpart.^^ Together with 
this observation, our result makes a three-dimensional 
cavity a promising component for hybrid-qubit architec- 
tures^ We note that Chiorescu et al. have made a 
similar observation in an electron spin ensemble utilizing 
a cavity ringing effect^ but did not address the A/TV- 
dependence of gc nor the parameter regime we report 
here. 

The microwave cavity we used is an X-band cylindrical 
dielectric ring resonator from Bruker (ER 4118X-MD5). 
Vc of this resonator is about 2 x 10 ~^ m^, giving an ex- 
pected gs/'^'T^ of 0.06 Hz. The refiection spectrum {Sn) 
from the resonator was measured with a network ana- 
lyzer (HP 8722D). As will be discussed later, in all the 
measurements, the sample position inside the cavity was 
carefully adjusted to maximize the coupling. 

Two materials were examined in this study; the first 
material is 2,2-diphenyl-l-picrylhydrazyl (DPPH) pow- 
der, a radical commonly used for field calibration of EPR 
spectrometers. It has an EPR linewidth of 2.1 G at room 
temperature. In this material, the spin concentration 
can be made high; our samples had a total number of 
spins A/'tot ranging from 1.0 x 10^^ to 2.4 x 10^^. As 
the spins are in thermal equilibrium at room tempera- 
ture, the number of polarized spins, which gives N in the 
present study, is given as A/'tot x (^ 9.8 GHz/2/cb 300 K), 
and varies from 7.8 x 10^"^ to 1.9 x 10^^. The second 
material is x-form lithium phthalocyanine (LiPc) pow- 
der electrochemically prepared from Li2Pc.^^ By reduc- 
ing the partial pressure of oxygen, this radical has a very 
narrow EPR linewidth, in our case around 83 mG at room 
temperature. On the other hand, N in these samples was 
much smaller than that of DPPH: less than 7.4 x 10^^ 
(including the polarization factor). The number of spins 
was determined with an X-band cw EPR spectrometer 
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FIG. 1. (Color online) (a) Reflection spectrum of DPPH with 
N = 7.8 X 10^^ in a logarithmic scale. The input microwave 
power was —20 dBm. The dashed line is a fit by Eq. ([Tj). 
The right panel shows the cross sections at 3470.9 G (black, 
on resonance) and 3466 G (gray), (b) Reflection spectrum of 
DPPH with AT = 1.9 X 10^^ The dashed lines are fits by 
Eq. Q, giving the splitting of 2^c/27r = 11.9 MHz. The cross 
sections at 3469 G (black, on resonance) and 3465 G (gray) 
are given in the right panel, (c) Change of n extracted from 
(a). The solid line is a Lorentzian fit. (d) gd"^^ of DPPH as 
a function of \fN . The triangular point (A) is obtained from 
the analysis in (c), whereas circular points (Q) are obtained 
from fits using Eq. Q. The point in gray corresponds to (b). 



(Bruker EMX), independently from the measurements 
described below. 

We first discuss the result from DPPH. Figure [2- (a) 
shows the reflection spectrum for the N = 7.8 x 10^^ 
sample. Here we observe features familiar in conven- 
tional EPR. As is well-known, the dispersion and absorp- 
tion components of the EPR signals are mathematically 
connected each other by the Kramers-Kronig relations.— 
The former appears as the shift of the cavity frequency 
across the resonance (at Bj. = 3470.9 G)^^ while the 
change in the cavity linewidth, or cavity quality factor 
Q, is a consequence of the latter. A model with two cou- 
pled oscillators describes the changes in cavity frequency 



00 and cavity half- width k as^ 
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where Kc is the cavity half- width far from the resonance, 
7s is the spin half- width, and A = mo{B — Br)/ h is 
the the field detuning. A fit to the cavity dip posi- 
tions in Fig. [I]- (a) gives ^c/27r = 1.15 MHz and 7s/27r 
= 2.85 MHz. Figure [I]-(c) plots as a function of ex- 
ternal magnetic fields. A fit using Eq. (j2j) gives ^c/27r = 
1.12 MHz, 7s/27r = 2.00 MHz, and Kc/^tt = 0.73 MHz. 
The values of gc obtained from the two fits agree well as 
expected. Those of 7^ slightly differ from each other but 
are consistent with that from cw EPR, in which the Gaus- 
sian full- width of 2.1 G corresponds 7s/27r = 3.5 MHz. 

As we increase the spectrum changes dramatically. 
When N is increased to four times that of the previous 
measurement, we start to observe two clearly separated 
dips, and the splitting becomes larger as we increase A^. 
Figure [B(b) shows the case for N = 1.9 x 10^^. The 
upper and lower branches uj± are fitted by the expression 
for the vacuum Rabi splitting 



, A , VA2+4gg 
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from which we deduce ^c/27r = 5.9 MHz. Figure [T]-(d) 
shows the V^- dependence of gc. The circular points 
represent samples in which the splitting was evident and 
from the fit we obtain ^s/27r = 0.043 Hz, which is in rea- 
sonable agreement with our estimation of 0.06 Hz. The 
triangular point is for the N = 7.8 x 10^^ sample that 
did not show splitting, but it is apparent that the value 
of gc extracted from the change in Q-factor also falls onto 
the fitted line to obey the \/]V- dependence. 

In theory, the observation of vacuum Rabi splitting 
and the realization of the strong coupling regime require 
the condition gc ^ i^c^ls or, similarly, the cooperativity 
C = gl/2K.cls > 1. This is barely satisfied in the A" = 1.2 
X 10^^ and 1.9 x 10^^ samples, where gc/27i = 4.9 MHz 
and 5.9 MHz are larger than both 7s/27r = 3.5 MHz and 
Kc/27r < 1 MHz (giving C ~ 4 and 7, respectively). On 
the other hand the N = 7.8 x 10^^ sample has gc larger 
than f<ic but smaller than 7^, and C of 0.2, consistent with 
the absence of the anticrossing. A question arises for the 
N = 3.1 X 10^^ sample, where gc/'^Tr = 2.0 MHz is com- 
parable with, but smaller than 7s/27r {C = 0.8) and yet 
we observed the splitting. This observation motivates us 
to explore a different parameter regime where gc is much 
smaller than either or 7^. This is possible with LiPc, 
in which, according to N available, ^c/27r will be less 
than 1 MHz while nd^ii can be larger than 1 MHz since 
this material is more lossy than DPPH. Owing to its nar- 
row EPR linewidth, 7s/27r is only 0.14 MHz, which may 
compensate the reductions in the other two parameters. 

Figure [2]- (a) demonstrates an anticrossing behavior in 
LiPc, where gd^i^ = 0.71 MHz and /^c/27r = 5.4 MHz 
are experimentally obtained. The appearance of two dips 
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FIG. 2. (Color online) (a) Reflection spectrum of LiPc with 
iV = 2.2 X 10^^. The cross sections are given for 3469.2 G 
(black, on resonance) and 3468.5 G (gray), (b) Simulation to 
(a). The same color bar as (a) applies. See the text for detail, 
(c) 5'c/27r of LiPc as a function of The square point (□) 
is omitted from the fit. (d) ^c/27r of the AT = 2.2 x 10^^ 
sample as a function of the sample position along the long 
axis of the cylindrical resonator. The origin of the horizontal 
axis is approximately the center of the resonator. The solid 
line is a fit with a sinusoidal function. In (c) and (d) the point 
in gray corresponds to (a). 



is obvious, although even at the resonance they are not 
fully separated (note that the data is given in dB units). 
As shown in Fig. [2]-(c), the a/]V- dependence of gc is con- 
firmed again with nearly the same gc as in the case of 
DPPH (0.043 Hz for DPPH and 0.049 Hz for LiPc). As 
for the = 7.4 x 10^^ sample, which was not used for 
the fitting, gc deviates from the V^-dependence. We at- 
tribute this as due to its very large sample volume com- 
pared with other samples, and the ac magnetic field the 
sample experiences is no longer uniform. We can esti- 
mate the distribution of the ac magnetic field by simply 
moving the sample across the cavity [Fig. [2]-(d)]. Indeed, 
it was necessary to check the position dependence in all 
samples to ensure that the coupling was maximized. Fit- 
ting the data with a sinusoidal function and assuming the 
same distribution in the N = 7.4 x 10^^ sample, we nu- 
merically average the coupling strength over the sample 
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FIG. 3. (Color online) (a) Simulation with gc/27v = 0.14 MHz 
(dc/js = !)• The dashed lines are Eq. Q. (b) gc/'^ir — 
0.085 MHz {gchs = 0.61). The dashed line is Eq. In 
both cases, all the parameters except for gc are the same as 
those for Fig. [1(b). 



length, and obtain gd'^^ = 1 MHz, which is reasonably 
close to the observed value of 0.9 MHz. 

We now discuss the condition to observe the splitting in 
more detail, using the following expression derived from 
a standard input-output formalism: -i^^ 



|5iii 



i(w - Wc) - Kc + gll {iA - 7^) 



(4) 



where Hie is the external loss. To demonstrate the validity 
of Eq. dll, we compute a spectrum using the parameter 
values gc/27r = 0.71 MHz, /^c/27r = 5.4 MHz, -fsf^n = 
0.14 MHz, and Kc = 0.99 x Kc [Fig. [2Kb)]. This repro- 
duces well the experimental spectrum in Fig. [2]- (a). Note 
that in this simulation the only adjustable parameter was 
and experimentally obtained values were used for gc, 
Kc, and 7s. Here, we see hZc > 9c > Js^ and C is only 0.3. 
The reason we still observe the splitting is understood 
through a straightforward algebraic analysis of Eq. (|4|). 
We differentiate it with respect to u to find local minima 
at A = 0. Assuming Hie/f^c ^ 1, we observe that two 
minima exist if — 7^(1 -\-4C){jg — 2^^) > is satisfied. 
For hZc large compared with gc and 7^, the condition is 
simplified to ^c/7s > Vv2 — 1 ^ 0.64. This indicates 
that the splitting may be observable even if gc is smaller 
than both 7^ and 

hhc 5 as confirmed by additional simula- 
tions shown in Fig. [3l As reducing only gc while all the 
other parameters fixed from Fig. [2]-(b), we observe that 
two branches merge into one at around gc/js = 0.64. We 
emphasize that this does not mean that the light-matter 
coherent information transfer is possible in this param- 
eter regime. For the coherent transfer to be possible, 
the condition gc ^ i^c^Js remains a compelling physical 
requirement. 

In conclusion, we have demonstrated the strong cou- 
pling between an electron spin ensemble and a three- 
dimensional cavity {gc > hCcHs) confirmed the >/iV- 
dependence of the collective coupling strength gc. For 
such a spin ensemble to be useful as a quantum memory, 
strong coupling and consequently entanglement with a 
single photon must be achieved, while the present exper- 
iments required many photons (10^^ or less) to spectro- 
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scopically resolve the anticrossing. In this respect DPPH 
suffers a relatively large inhomogeneous linewidth, which 
results in a short storage time unless multiple microwave 
pulses for decoupling/refocusing are applied. We have 
also observed an anticrossing behavior in the spectrum 
with the cooperativity C smaller than 1, and derived 
the condition for observing the splitting. An advan- 
tage of our system over a cavity QED counterpart lies in 
the simplicity in the experimental setup; we may readily 
choose other spin systems with different EPR linewidths 
(to control 75) and spin densities (to control A/"), and N 
may also be controllable by the temperature or by means 
of dynamic polarization schemes. Therefore, we believe 
that our observation will broaden the applicability of the 
physics of collective light-matter coupling to wider pa- 
rameter regimes. 
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